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Dynamical model of electroweak pion production in the resonance region 
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In this report, we will briefly review the dynamical model of pion electroweak production reactions in the A 
resonance region and report on our study of neutrino-nucleus reactions based on this model. 



1. Introduction 

It is well recognized that the precise knowledge 
of neutrino-nucleus reactions is crucial in analyz- 
ing neutrino oscillation experiments pQ . For neu- 
trino reactions around 1 GeV, quasi-elastic scat- 
tering and pion production processes are the main 
reaction mechanisms. The pion production takes 
place through both non-resonant and resonance 
processes, and in the latter the A 33 resonance 
plays a central role for neutrinos in the 1 GeV 
region. Since experimental data on the neutrino- 
nucleon reactions in the resonance region are not 
as extensive as electron scattering data, a theoret- 
ical study of the weak pion production amplitude 
will be valuable for the study of neutrino-nucleus 
reactions |2I3I4I5I6I7| . 

In the recent years, extensive studies on the 
electron- and photon-induced meson production 
reactions in the GeV region are in progress to 
investigate the nucleon resonance properties [HJ- 
The objective of the N* study is to understand 
the non-perturbative dynamics of QCD by test- 
ing the resonance properties against the predic- 
tion of QCD-inspired models and/or lattice sim- 
ulations. Two of the present authors with their 
collaborators have developed a dynamical model 
for describing photo- and electro-production of 
pions off the nucleon around the A resonance re- 
gion \21A\ . It has been shown that including the 
pion cloud effect in theoretical analyses can re- 
solve a long-standing puzzle that the JVA mag- 
netic dipole transition form factor Gm predicted 



by the quark model is about 40% lower than the 
empirical value. Furthermore, the predicted E2 
Ge and C2 Gc form factors show pronounced Q 2 
dependence due to the pion cloud effects, which 
suggests the deformation effects in the NA tran- 
sition. This model is further extended to investi- 
gate neutrino reactions, which involve the axial- 
vector responses of hadrons jlloj . The NA axial 
vector form factor Ga is found to contain large 
meson cloud effects. 

The dynamical model starts from the non- 
resonant meson-baryon interaction and the res- 
onance interaction, and the unitary amplitudes 
are obtained from the scattering equation. Fairly 
consistent descriptions of all the available data in 
the A resonance region have been obtained from 
the dynamical model. The well tested dynamical 
model for the single-nucleon case will be a good 
starting point for investigating the reaction dy- 
namics of neutrino- nucleus reactions. 

In this paper we give a short review of our dy- 
namical model for meson production reactions in 
the delta resonance region. Subsequently we show 
the results of our first application of this model 
to neutrino-nucleus reactions in the delta produc- 
tion region. 

2. Dynamical model 

We briefly describe our dynamical model of 
pion photoproduction. We start from the inter- 
action Hamiltonian (H) for the N,A,ir,p and u> 
fields. By assuming that only 'few-body' states 
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are active in the energy region we are interested 
in, we can simplify the intrinsic many-body prob- 
lem. The effects due to many-body states are ab- 
sorbed in effective interaction operators. In our 
work, the effective 'few-body' Hamiltonian -ff e // 
is derived from H using the unitary transforma- 
tion method (Hj- The resulting effective Hamil- 
tonian H e ff is defined in a subspace spanned by 
irN, A, 7 AT states and has the following form: 

Heff = Ho + TttAt^a + T 7 at + _ > a + VnN + %ZV- (1) 

Here v^n and v 7 jv are non- resonant interactions, 
which consist of the u- and s-channel nucleon ex- 
change and t-channel pion and vector meson ex- 
change processes shown in Fig. 1(a). The exci- 
tation of A is described by r T ^„A and T 7 jv^a 
shown in Fig. 1(b). An important feature of our 
approach is that the effective Hamiltonian is en- 
ergy independent. Hence the unitarity of the re- 
sulting amplitude is automatically satisfied. Fur- 
thermore, the non-resonant interactions v^n and 
i> 7 jv are derived from the same unitary transfor- 
mation and hence 7riV and 7-/V reactions can be 
described consistently. 



(a) 



(b) 



Figure 1. Graphical representation of jN 
and 7 TV — * ttN interactions. 



From the effective Hamiltonian, it is straight- 
forward to derive a set of coupled equations for 
irN and 7 TV reactions. The resulting pion photo- 
production amplitude is given as 



T 77r = <7riV|e- J em |iV> 



^on+ ^rgyrsaf . (2) 
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where W is the invariant mass of the irN sys- 
tem. The first term is the non-resonant amplitude 



i 77r calculated from the non-resonant interactions 
u^jv and i> 7 at. The second term in Eq. © is the 
resonant amplitude. The self-energy E is defined 
as 

S(W0 = T KN ^ A G t A ^ N {W). (3) 
The 'dressed' 7-/V — > A vertex is defined by 

f jN ^ A {W)=T jN ^ A + J T nN ^ A G t 77T (W). (4) 

An important consequence of the dynamical 
model is that the influence of the non-resonant 
mechanisms on the resonance properties can be 
identified and calculated explicitly. The 'bare' 
resonance vertex T 7 n^a is modified by the non- 
resonant meson cloud i 77r to give the dressed ver- 
tex F 7 AT^A- 

3. Pion electroweak production and NA 
form factors 

The dynamical approach described in the pre- 
vious section was used to study n-N scatter- 
ing, pion photoproduction and pion electropro- 
duction reactions |2I3| . The model was subse- 
quently extended to investigate neutrino induced 
reactions | II") j . Since the weak currents in the 
standard model are closely related to the electro- 
magnetic current, it is straightforward to extend 
the dynamical model of pion electroproduction to 
weak pion production reactions. The electromag- 
netic current (jff M ), the weak charged current 
(j^ c ) and the weak neutral current (jf? ) can 
be expressed in terms of the vector current 



and the axial-vector current as 
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(5) 
(6) 
(7) 



Using CVC, V^- 2 can be obtained from V$ by 
isospin rotation. Guided by the effective chiral 
Lagrangian method and using the unitary trans- 
formation method, we can construct A^ for pion 
production. The resulting A^ consists of the 
nucleon-Born term, rho-exchange and delta ex- 
citation terms. The 'bare' NA magnetic and 
electric form factors Gm(0) and Ge{0) are ob- 
tained from the analysis of the (7, it) reaction. 
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We assume the Siegert theorem to determine the 
Coulomb quadrupole form factor G*c(0) from the 
transverse electric form factor Ge(0). The Q 2 
dependence of the form factors are assumed to be 

G a (Q 2 ) = G a (0)R SL (Q 2 )G D (Q 2 ), (8) 

where Gd = 1/(1 + Q 2 /m 2 /) 2 is the dipole form 
factor of the proton with m 2 , = 0.71 GeV 2 . The 
correction factor R$l is defined as 

Rsl{Q 2 ) = (l + aQ 2 )exp(-6Q 2 ), (9) 

where a = 0.154 GeV -2 and b = 0.166 GeV -2 
have been deduced by fitting the (e, e'ir) data 
at Q 2 = 2.8 and 4 GeV 2 10 . For the iVA ax- 
ial vector form factor Ga(Q 2 ), we also assume 
the form given in Eq. © with G D {Q 2 ) replaced 
by the nucleon axial-vector dipole form factor 
with — 1-02 GeV. We further assume the 
SU(6) quark model relation between Ga for the 
NA transition and gA of the nucleon. Having 
determined the axial NA form factor, we have 
no further adjustable parameter in our model to 
describe neutrino-induced pion production reac- 
tions. 

In Fig. 2, the energy dependences of the 
p(e, e'n + )n and p{v, e~n + )p cross sections are 
shown at Q 2 — 0.1(GeV/c) 2 and Ei epton = 1 
GeV. The solid lines show the full results and 
the dashed lines show the contribution of the res- 
onance amplitude which is the second term of 
Eq. 2. Clearly the non-resonant amplitude plays 
a crucial role. It is well known that the Kroll- 
Ruderman term dominates low-energy charged- 
pion production reactions. We note that even in 
the p{v, e~ir + )p reaction, where only the isospin 
3/2 amplitude contributes, the non-resonant con- 
tribution to the cross section amounts to 15% at 
W = 1.2 GeV. Thus the appropriate treatment of 
the non-resonant mechanisms plays an important 
role in interpreting the NA transition form fac- 
tors G a (Q 2 ) and also in determining the neutrino- 
nucleon reaction amplitudes that are needed for 
calculating the neutrino-nucleus reaction ampli- 
tudes. 

The pion production mechanism of our dynam- 
ical model may be tested by comparing the the- 
oretical and experimental values of the cross sec- 
tions for a wide energy ranges also useful. For 
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Figure 2. Cross section da/dW/dQ 2 of 
p(e,e'ir + )n (left) and p(v, e~ir + )p (right) reac- 
tions. 



instance, the longitudinal-transverse interference 
term ctlt' for a polarized electron is given by in- 
terference between the real and imaginary parts 
of the amplitude and hence it is sensitive to 
the non-resonant mechanism. As an example, 
we compare in Fig. 3 our predictions for the 
e+p — > e'+p+n reaction with the Jlab data [11| . 
The cross section for virtual-photon pion produc- 
tion can be expressed as 
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(10) 



The transverse cross section <tt is mainly deter- 
mined by the contribution of the magnetic dipole 
form factor Gm ■ Interference between the trans- 
verse and longitudinal currents (<jlt) is sensitive 
to the product GcGm- Our model gives a reason- 
able description of the angular distribution from 
threshold to the delta resonance region. Reason- 
able agreement has also been seen in the compar- 
ison of our dynamical model predictions with the 
extensive pion production data in the delta reso- 
nance region obtained at LEGS, Mainz, Jlab and 
MIT-Bates. 

We compare in Fig. 4 the calculated Q 2 de- 
pendence of the cross section for the + p — > 
fi~ + 7r + + p reaction with the ANL ^21 an d 
BNL ^2] data. Here we take into account the 
variation of the neutrino flux in the experiment 
and the finite mass of the muon. The results 
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Figure 3. p(e,e'TT°)p cross section at Q 2 
0.4(GeV/c) 2 . 
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Figure 4. Differential cross section da/dQ 2 of p + 
v u — * M _ + 7r+ + V reaction 



and the dressed axial- vector form factor, Ga, ex- 
tracted from our present model. The pion cloud 
effect is seen to be essential in explaining the em- 
pirical values of Gm directly extracted from the 
data. The pion cloud effects are sizable also for 
Ga- The dynamical model explains why most of 
the quark models underestimate Gm and Ga- 



of our full calculation (solid line) agree reason- 
ably well with the data both in magnitude and 
Q 2 dependence. The contribution of the axial 
vector current (dashed curve) and the vector cur- 
rent (dot-dashed curve) have rather different Q 2 
dependences in the low Q 2 region and interfere 
constructively with each other. Since the vector 
current contributions are highly constrained by 
the (e, e'n) data, the results in Fig. 4 suggest 
that our model for the axial vector current reac- 
tion is consistent with the data. This in turn im- 
plies that the axial-vector NA coupling constant 
predicted by the SU(6) quark model is consistent 
with the existing data of neutrino-induced pion 
production reactions in the A resonance region. 
A precise measurement of parity- violating asym- 
metry of inclusive p(e, e') can be used to improve 
the determination of the axial NA form factor |S] . 

Finally, we explore the effects of the pion cloud 
on the NA transition form factors. Fig. 5 displays 
the dressed magnetic dipole form factor, Gm, 
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Figure 5. GM(left) and G^right). 



4. Neutrino-nucleus reaction 

We now apply the dynamical model explained 
in the previous section to neutrino-nucleus reac- 
tions in the GeV region. The charged-current 
neutrino-nucleus reaction cross section can be 
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written as 

da p[ G 2 F cos 2 8 C 

dE'dfl' ~ 



Pi 



8tt 2 



(11) 



where L^ v is the lepton tensor. The hadron ten- 
sor W^ v is given as 



f) 



<f\j»\i><f\r\i>*, 



(12) 



where Pij is the initial or final hadron momen- 
tum, q a momentum transfer. 

As a first step of our neutrino-nucleus reac- 
tion study, we consider the following points: (1) 
Fermi-averaging of the elementary amplitude is 
taken into account by calculating the reaction en- 
ergy W = \/(q + Pn) 2 of the elementary process 
from the nucleon momentum p^ inside nucleus 
and momentum transfer q. (2) We also take into 
account Pauli effects on the final nucleon after 
pion production. We require that the momentum 
of the final nucleon = \p^ + q — k\ be larger 
than the Fermi momentum. This correction is 
calculated for each pion momentum k and initial 
nucleon momentum p^. Fig. 6 shows the cross 



section for the 



+ TT + N+(A = 11) 



reaction divided by mass number for E v — 1 GeV 
and for the two values of the lepton scattering 
angle, 10° and 30°. The dash-dotted curve shows 
the average of the cross sections for the proton 
and neutron targets, (a p +a n )/2. Fermi averaging 
is included in the dotted curve, and furthermore 
the Pauli-effect is included in the solid curve. As 
discussed in Ref. [Hj, the Pauli effect is apprecia- 
ble at low momentum transfers, i.e. at forward 
lepton angles. In evaluating the Pauli effect, we 
take full account of the pion angular distribution 
as predicted by the dynamical model. 

To examine to what extent our treatment can 
describe the data, we have studied inclusive elec- 
tron scattering on 12 C. The quasi-elastic pro- 
cess is evaluated using the relevant formula in 
Ref. 14 . The pion production process is calcu- 
lated as described above. Here the nuclear corre- 
lation is taken into account by using the structure 
function given in Ref. 15 . In Fig. 7, our results 
are compared with the data for E e = 0.96 and 




Figure 6. Cross section of v + 12 C — > e +tt + X. 



1.1 GeV at 9 = 37.5° The magnitude of the 
calculated pion production cross section is com- 
parable with the data, while our result (solid line) 
tends to overestimate the quasi-free cross section. 
Furthermore the 'dip' region between the quasi- 
free and delta-excitation peaks is not explained 
within the model. 

Recently there have been detailed studies of 
nuclear effects in inclusive neutrino-nucleus reac- 
tions |17I18| . Meanwhile, in the present first at- 
tempt to apply the dynamical model of Refs. 
I3I4I5| to neutrino-nucleus reactions, we have in- 
cluded nuclear-medium effects only by incorpo- 
rating the Fermi-averaging and Pauli effects in 
the way described above and by introducing the 
structure function. It is obviously important to 
investigate any possible nuclear effects beyond the 
level considered in our present work. 
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Figure 7. Inclusive cross section of 12 C + e — * 
e' + X at E e = 0.96 GeV (left) and E e = 1.1 GeV 
(right). 
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5. Summary 

We have developed a dynamical model of elec- 
troweak pion production reactions. Most of the 
available data of pion electroproduction and neu- 
trino reaction on the nucleon in the A region 
can be described reasonably well by the model. 
The iVA transition form factors extracted from 
the data exhibit a large contribution of the pion 
cloud effects to the NA form factors, and these 
effects explain why the simple quark models fail 
to explain the magnitude and the Q 2 dependence 
of the empirical form factors. We have studied 
neutrino-nucleus reactions using the dynamical 
model that has been well tested by the electron 
scattering data. We have examined the effects of 
Fermi motion, Pauli blocking and nuclear correla- 
tions using the dynamical model. It is shown that 
these effects are important and significantly im- 
prove the description of the electron-nucleus scat- 
tering around lGeV. However further studies on 
the propagation of the A-particle in nuclei seem 
to be needed. 

Acknowledgments 

This work is supported by the Japan Society for 
the Promotion of Science, Grant-in- Aid for Scien- 
tific Research(c) 15540275, by the U.S. National 
Science Foundation, Grant Nos. PHY-0140214 
and PHY-0457014, and by the U.S. Department 
of Energy, Nuclear Physics Division Contract No. 
W-31-109-ENG-38. 



6. E. A. Paschos, Ji- Young Yu and M. Sakuda, 
Phys. Rev. D 69 (2004) 014013. 

7. O. Lalakulich and E. A. Paschos, Phys. Rev. 
D 71 (2005) 074003. 

8. V. D. Burkert and T. -S. H. Lee, Int. J. Mod. 
Phys. E 13 (2004) 1035. 

9. M. Kobayashi, T. Sato and H. Ohtsubo, Prog. 
Theor. Phys. 98 (1997) 927. 

10. V. V. Frolov et al., Phys. Rev. Lett 82 (1999) 
45. 

11. K. Joo et al., Phys. Rev. Lett. 88 (2002) 
12001. 

12. S. J. Barish et al., Phys. Rev. D 19 (1979) 
2521. 

13. T. Kitagaki et al., Phys. Rev. D 42 (1990) 
1331. 

14. R. A. Smith and E. J. Moniz, Nucl. Phys. B 
43 (1972) 605. 

15. O. Benhar, A. Fabrocini, S. Fantoni and I. 
Sick, Nucl. Phys. A 579 (1994) 493. 

16. R. M. Sealock et al., Phys. Rev. Lett. 62 
(1989) 1350. 

17. O. Benhar et al., Phys. Rev. D 72 (2005) 
053005. 



18. J. Nieves et al, nucl-th/0503023 



REFERENCES 

1. J. G. Morfin and M. Sakuda and Y. Suzuki 
(eds), Proceedings of the First International 
Workshop on Neutrino-Nucleus Interactions 
in the Few GeV Region (NuIntOl), Nucl. 
Phys. B(Proc. Suppl.) 112 (2002). 

2. T. Sato and T.-S. H. Lee, Phys. Rev. C 54 
(1996) 2660. 

3. T. Sato and T.-S. H. Lee, Phys. Rev. C 63 
(2001) 055201. 

4. T. Sato, D. Uno and T.-S. H. Lee, Phys. Rev. 
C 67 (2003) 065201. 

5. K. Matsui, T. Sato and T. -S. H. Lee, Phys. 
Rev. C 72 (2005) 025204. 



